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ABSTRACT 

We have obtained deep optical spectrophotometry of 16 planetary nebulae in M33, 
mostly located in the central two kpc of the galaxy, with the Subaru and Keck tele- 
scopes. We have derived electron temperatures and chemical abundances from the 
detection of the [O ill] A4363 line for the whole sample. We have found one object with 
an extreme nitrogen abundance, 12+log(N/H) = 9.20, accompanied by a large helium 
content. After combining our data with those available in the literature for PNe and 
H II regions, we have examined the behavior of nitrogen, neon, oxygen and argon in 
relation to each other, and as a function of galactocentric distance. We confirm the 
good correlation between Nc/H and 0/H for PNe in M33. Ar/H is also found to cor- 
relate with 0/H. This strengthens the idea that at the mctallicity of the bright PNe 
analyzed in M33, which is similar to that found in the LMC, these elements have not 
been significantly modified during the dredge-up processes that take place during the 
AGB phase of their progenitor stars. We find no significant oxygen abundance offset 
between PNe and Hii regions at any given galactocentric distance, despite the fact 
that these objects represent different age groups in the evolution of the galaxy. Com- 
bining the results from PNe and Hii regions, we obtain a representative slope of the 
ISM a-element (O, Ar, Ne) abundance gradient in M33 of -0.025 ± 0.006 dexkpc"^ 
Both PNe and H ii regions display a large abundance dispersion at any given distance 
from the galactic center. We find that the N/0 ratio in PNe is enhanced, relative to 
the Hii regions, by approximately 0.8 dex. 

Key words: galaxies: abundances - galaxies: ISM - planetary nebulae - galaxies: 
individual: M33. 



1 INTRODUCTION 

Planetary nebulae (PNe) can be used as test particles in 
the study of the chemical enrichment of galaxies and as fos- 
sil records of the nucleosynthesis that took place in previ- 
ous generations of stars. While information on present-day 
chemical abundances of star-forming galaxies are gathered 
from the spectra of H ll regions and massive stars, PNe pro- 
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vide the chemical makeup of older progenitor stars. More- 
over, imprinted on their emission-line spectra are the sig- 
natures of chemical elements synthesized during the AGB 
phase of low- and intermediate-mass stars. Various processes 
occurring in these stars, such as third dredge-up and hot bot- 
tom burning, can strongly modify the helium, nitrogen and 
carbon abundances of the outer layers, which are eventually 
blown out to form the observed PNe. The effects on other 
elements, such as oxygen and neon, can be more subtle, and 
are currently the su bject of lively debate (|Stasinskal l2008l : 
ICristallo et aLll2009^ . 

PNe represent an alternative to H ll regions in the mea- 
surement of galactic abundance gradients, allowing tests of 
the evolution predicted by models of gal axy assembly (e.g. in 
the c ontext of an inside-out formation, iPrantzos fc Boissiej 
l200d ). In addition, PNe can be used to test H ll region abun- 
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dances in cases where the latter are expected by theoreti- 
cal calculations to be affected by systematic errors. These 
tests are particularly important in the central, metal-rich re- 
gions of spiral galaxies, where the high metal content could 
lead to systemat ic underestimates of the nebular metallicity 
(|Stasihskal |2005| ) . This originally motivated our new inves- 
tigation of a sample of luminous PNe in the central region 
of the galaxy M33. Despite the faintness of the targets in 
this external galaxy, compared to objects in the Milky Way, 
the uncertainty on the relative galactocentric distances, that 
afflicts studies in the Galaxy, is virtually removed. This of- 
fers a clear advantage in the determination of abundance 
gradients, and consequently in deriving information on the 
chemical evolution of the host galaxy, as well as provid- 
ing empirical constraints on the metallicity dependence of 
extragalactic distance indicators (e.g., the Cophoi d period- 
lumin osity relation: iKcn nicutt ct al. 1998, Scowc roft et al.l 
I2OO9D . 

In this paper we present new spectroscopic observa- 
tions of 16 PNe in the neighboring galaxy M33, obtained 
with the aim of measuring chemical abundances from the 
[O III] A4363 diagnostic line in the central few kpc of this 
galaxy. This classical 'direct' method is reckoned to yield 
the most accurate nebular abundances, unaffected by sys- 
tematic uncertainties and/or strong-line method calibration 
issues, although in extragalactic PN and H 11 region work its 
application is made difficult, especially at high metallicity 
or large redshift, by the weakness of the A4363 auroral line. 
Combining our new data with existing observations in the 
literature we discuss the abundances of He, O, Ne, N and 
Ar in relation to each other, and as a function of galacto- 
centric distance, in order to provide a clearer picture of the 
chemical composition and evolution of M33. The number of 
recent studies o n the present-d ay chemi cal composition of 
ioniz e d nebulae dCrockett et al. 2006; Rosolowskv &i Simonl 
l2008l : iMagrini StangheUini. fc ViUaverl 120091 ') and massive 
stars ((u et al.ll2009l T in M33 attests to the key role played by 



Table 1. Log of the observations (Subaru/FOCAS). 



this Milky Way neighbour in the investigation of the chem- 
ical evolution of spiral galaxies. 



2 OBSERVATIONS AND DATA REDUCTION 

Deep multi-slit spectra of PNe in the central region of 
M33 were obtained at the 8.2m Subaru Telescope equipped 
with the Faint Object Camera and Spectrograph (FOCAS, 
iKashikawa et al" I I2OO2I ) on 2007 October 9-10. Narrow-band 
[O III] A5007 images obtained with the same in strument were 
used, toget her with the PNe coor dinates from lMagrini et al] 
l|200lh and lCiarduUo et all (|2004 ). to select targets and pre- 
pare the multi-object masks in two circular fields of 3' radius. 
One of the fields contained the nucleus of the galaxy (the 
field center is approximately l'.5 arcmin SE of the galaxy 
center), while the second one was centered 7' to the N of 
the nucleus (coordinates of the field centers are given in Ta- 
ble [l]). 

For the spectroscopic observations we used 1'.' 2- wide 
slits and a combination of three different grisms (300R in 
second order, 300B and VPH650) to secure a continuous 
wavelength coverage from the near-UV (in most cases in- 
cluding the [0 11] A3727 line) to the red (up to [Arm] A7135), 
with a maximum spectral resolution of 4.5 A (300R grism in 



Field 


R.A. 

Dec. 

(J2000) 


Exposure 
time 

(s) 


Grism 


1 


0l''33™54!14 


3 X 1800 


300R 




30°38'16'.'5 


1 X 1200 


300B 






2 X 1200 


VPH650 


2 


0l''34™03!61 


3 X 1800 


300R 




30°45'43('6 


1 X 900 


300B 






3 X 1200 


VPH650 



the blue) and 4.0 A (VPH650 grism in the red). Seeing con- 
ditions during the exposures presented here varied between 
0.8 and 1.0 arcsec during the first night, and between 0.5 
and 0.7 arcsec during the second night. Additional spectra 
were obtained in 2008 under worse atmospheric conditions 
(seeing up to 2 arcsec), but due to their poor quality they 
were not included in the analysis. The spectrophotometric 
standards G 191-B2B, GD 71 and LDS 749B were observed 
to flux calibrate the PNe spectra. Table [1] summarizes the 
exposure times used with the various grisms. 

The data reduction was carried out by means of stan- 
dard iRAlj^ tasks, and included bias and flat-field cor- 
rections, wavelength calibration from thorium-argon lamp 
expo sures, cosmic-ray r emoval with the L. A. COSMIC rou- 
tine (|van DokkumllioOll ). flux calibration and spectra co- 
addition. The subtraction of the background was carried 
out locally for each individual PN, from the flux measured 
within the corresponding slitlet, each of which was suffi- 
ciently long (~20 arcsec) to allow a good sampling of the 
background (a linear fit to the background measured in the 
vicinity of each PN was sufficient). This step is crucial for 
the correct measurement of the low-excitation lines, which 
can be present in emission also in the diffuse ISM. It is worth 
pointing out in this context that in general the sky subtrac- 
tion obtained from slit spectra is known to be superior to 
what can be achieved with fiber-fed spectrographs, being 
less prone to systematic errors. 

PNe spectra were also obtained on 2004 September 22- 
23 and October 21 with the Keck I telesc ope and the Low - 
Resolution Imaging Spectrometer (LRIS, lOke et al] Il995l ) 
as part of the spectrosc opic survey of Hll re g ions in the 
southern haff of M33 bv iRosolowskv fc Simonl (120081 ). The 
Keck sample extends the galactocentric distance coverage 
of our sample to larger radii. The blue spectra presented 
here, extending from [0 11] A3727 to [O ill] A5007, were ob- 
tained with the 600/4000 grism, providing '-^5A resolution 
with the l"-wide slits used in the multi-object masks. More 
details on the reduction of t he Keck dataset can be found in 
IRosolowskv fc Simonl (|2008l . =RS08). 

For our chemical abundance analysis we retained only 
spectra for which we had a reliable detection of the Te- 
sensitive line [O ill] A4363. This left us with 10 and 8 PNe 
from the Subaru and the Keck observations, respectively, 



^ IRAF is distributed by the National Optical Astronomy Observa- 
tory, which is operated by the Association of Universities for Re- 
search in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Spectra of PN 18, PN 55, PN 67 and PN 68. Zoomcd-in views are provided in the insets to better display the weakest features. 
The main emission hnes are identified. 



with two targets in common between the two datasets. Only 
two Subaru targets were excluded from the analysis, because 
their [O ill] A4363 line could not be measured. The Subaru 
slit masks also included a few H ll regions. For three of them 
(B 90 in Field 1, B 72 and B 302 in Field 2; we ado pt the H ll 
region identification from I Boniest eix et al.) 1 19741 ) we had a 
good quality measurement of the electron temperature from 
the [O III] A4363 auroral lines. For the comparison between 
PNe and Hll region chemical abundances (Section 5) we 
retain these additional objects. For B 72 and B 302, to our 
knowledge, no previous direct abundance determinations ex- 
ist in the literature. 

Celestial coordinates and galactocentric distances for 
the final sample of PNe and H ll regions are summarized in 
Table [2] where the PNe coordinates and identifications are 



taken from the catalog by ICiarduUo et al] (|2004l ). For the 
calculation of the de-projected galactocentric distances we 
adopted a disk inclinat ion angle i — 56° and a position an- 
gle e = 23°, as used bv lCiarduUo et all (|2004l '). The position 
of the galaxy nucleus was s et at RA(J2000) = 0l''33'"50:92, 
DEC(J2000) =30°39'36'.'8 (|Massev et al.lll996l ). For consis- 
tency with previous nebular abundance i nvestigations, we 
adopt ed the Cepheid distance of 840 kpc bv lFreedman et aU 
(|200lD . even though more recent st udies suggest consid- 
erably larger values (up to 968 kpc: iBonanos et al. I l2006l . 
lU et al.llioogi '). In Fig. [TH we show the spectra of the 16 
PNe included in this work. It can be seen that none of the 
targets has detectable WR-like features. 

As Table[2]shows, our sample includes 5 PNe with galac- 
tocentric distance i? < 1 kpc, and 11 are within 2 kpc from 
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Figure 2. Spectra of PN 73, PN 76, PN 78 and PN 83. Zoomcd-in views are provided in the insets to better display the weakest features. 
The main emission lines are identified. 



the galactic center (the 3 Hll regions observed by us also 
lie in the central 2 kpc). This central region of the galaxy 
is under-represented in emission-line studies of the chemical 
composition of M33, because of the observational difficul- 
ties related to the measurement of the faint auroral lines 
in metal-rich environments. We also point out that the vast 
majority of the targets are expected to be part of the disk of 
M33. From the measu rement of the kinematics of 140 PNe 
ICiarduUo et al.l ()2004l ) determined that only two might be- 
long to the halo of the galaxy, rather than the disk. One of 
them, PN67, is included in our sample. Lastly, it is worth 
pointing out that the PNe observed as part of our program 
occupy the brightest 1.5 msooT magnitudes in the luminosity 
distribution of M33 planetaries. 

The measurement of the emission line strengths was 



carried out with the splot program in iraf. The extinc- 
tion c(H/3) was derived iteratively from the Balmer decre- 
ment , assuming case B Hi line ratios jHummer fc StorevI 
Il987h valid for the electron temperatures determined in Sec- 
tion [3] The reddening-corrected line fluxes, normal ized to 
I(H/3) = 100, and calculated using the ISeatonI (|l979l ') inter- 
stellar reddening law, are presented in Tables [3] and 2] The 
extinction coefficient c(H/3) appears in the last column of 
Table O The line flux errors account for the uncertainties in 
the placement of the continuum, in the flat fielding and the 
flux calibration, and in the extinction coefficient. 

In two cases (PN55 and 111) our spectra did not cover 
the bluest wavelengths, hence the entries for [O ll] A3727 
and/or [Nein]A3868 are missing in Table [3] Similarly, the 
Keck spectra did not extend beyond [O ill] A5007, there- 
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Figure 3. Spectra of PN 86, PN 90, PN 94 and PN 99. Zoomcd-in views are provided in the insets to better display the weakest features. 
The main emission hnes are identified. 



fore all the line fluxes at longer wavelengths are missing 
for the corresponding targets in Table [l] In these tables 
we report the line fluxes measured independently from the 
Subaru/FOCAS and Keck/LRIS spectra for the two tar- 
gets (PN 76 and PN 78) in common. This offers us a way to 
check the consistency of the flux intensities between the two 
samples, obtained from different instruments and reduced 
independently from each other. Considering the l-cr errors, 
we flnd quite a good agreement between the independent 
measurements. 

We have com pared our fluxes with those of 
iMagrini et al] (|2009l . =M09) for the 11 PNe in com- 
mon. The comparison is shown in Fig. [S] and it includes the 
O, Ne, N, S and Ar metal lines, as well as Hel AA4471, 5876 
and HellA4686. We show in colour the points representing 



O lines ([Oii]A3727, [Oiii]A4363 and [O iii] A5007), used 
to derive the O abundance. The agreement is quite good 
for strong lines, and worsens with decreasing line strength. 
While most of the points lie within the 0.2 dex scatter 
region from the line of equality, there are several emission 
features that display a larger discrepancy. This is especially 
the case of the [S ll] AA6717, 6731 lines, which are often 
quite faint in the spectra of PNe at the distance of M33, 
and which appear to be systematically too strong in the 
M09 sample compared to ours (open red squares in Fig. [SJ. 
A possible explanation is that the M09 spect ra are affected 
by w eak emission from the diffuse ionized gas (|Haffner et al.l 
120091 ). which we clearly detect in our 2-D spectra and which 
might be not properly subtracted in the work by M09, who 
used fibers centered in low emission regions of the galaxy 
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Figure 4. Spectra of PN 105, PN 108, PNlll and PN 113. Zoomcd-in views arc provided in the insets to better display tlie weakest 
features. The main emission lines are identified.. 



for the sky subtraction. The diffuse emission would be most 
important for faint, low-excitation metal lines, such as those 
of [Sll], and the effects of an improper subtraction on the 
PN spectra would increase with decreasing line strength, as 
observed in Fig. [S] 

We also found a large systematic difference concerning 
the extinction coefficient, our values being lower than those 
reported by M09 by an average factor of 2.5. We point out 
that our Balmer line intensity ratios agree with the theo- 
retical values, while in several cases the H7/H/3 and H5/H/? 
ratios determined by M09 are significantly below the case B 
values. The means we derive from their data are 0.403±0.102 
and 0.220 ± 0.057, respectively (case B values are 0.469 and 
0.259 for Te = 10* K, Ne = 10^ cm"^). This suggests a sys- 
tematic over-estimate of the extinction coefficient. 



The mean reddening for our PNe sample corresponds 
to E{B - V) = 0.19 ± 0.03. This is similar to other spectro- 
scopic determinations of the reddening from studies of the 
young populations of M33, including H 1 1 regions (e.g. RS08 : 
E{B-V) = 0.20±0.03) and young stars jMassev et al.lll995l : 
eIb -V) = 0.13 ±0.01). The possibility that the PNe have 
reddening values that are much larger than those measured 
in young star-forming regions seems unlikely, even though a 
somewhat higher local reddening could be expected for PNe, 
since the progenitor AGB stars can produce dust shells. 



3 CHEMICAL ABUNDANCES 

For the derivation of the ionic abundances it is necessary 
to obtain the electron temperature Te of the emitting gas. 
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Table 2. PNe and Hll region sample. 



Number 10+ R.A. (J2000)t Dec. (J2000)t B} Telescope* c(H/3) 

h m s oil, (j^p^-) 













PNe 








1 


PN 18 


01 


33 


06.11 


30 31 04.5 


3.73 


K 


0.33 ±0.08 


2 


PN55 


01 


33 


40.16 


30 37 49.5 


0.88 


K 


0.12 ±0.09 


3 


PN67 


01 


33 


48.27 


30 33 15.7 


1.71 


K 


0.37 ± 0.09 


4 


PN 68 


01 


33 


48.57 


30 35 47.9 


1.01 


S 


0.50 ±0.13 


5 


PN 73 


01 


33 


51.06 


30 45 38.8 


1.70 


S 


0.00 ±0.09 


6 


PN 76 


01 


33 


52.77 


30 37 38.9 


0.64 


S,K 


0.03 ±0.08 


7 


PN 78 


01 


33 


54.68 


30 36 05.7 


1.17 


S,K 


0.27 ±0.14 


8 


PN 83 


01 


33 


57.18 


30 36 47.6 


1.14 


S 


0.02 ±0.11 


9 


PN 86 


01 


33 


59.97 


30 40 28.7 


0.74 


s 


0.31 ±0.10 


10 


PN 90 


01 


34 


02.57 


30 40 00.5 


1.00 


s 


0.16 ±0.09 


11 


PN 94 


01 


34 


03.54 


30 39 15.6 


1.17 


s 


0.05 ± 0.09 


12 


PN 99 


01 


34 


06.56 


30 48 23.3 


2.30 


s 


0.47 ±0.13 


13 


PN105 


01 


34 


11.91 


30 45 00.9 


1.87 


s 


0.00 ± 0.09 


14 


PN108 


01 


34 


13.40 


30 33 50.7 


3.05 


K 


0.15 ±0.12 


15 


PNlll 


01 


34 


14.79 


30 31 49.6 


3.62 


K 


0.22 ±0.10 


16 


PN113 


01 


34 


15.48 


30 32 20.3 


3.55 


K 


0.28 ± 0.08 












H 11 regions 








17 


B90 


01 


34 


04.34 


30 38 12.9 


1.40 


s 


0.22 ±0.07 


18 


B302 


01 


34 


06.92 


30 47 26.1 


2.09 


s 


0.25 ±0.08 


19 


B72 


01 


34 


01.21 


30 43 56.8 


1.20 


s 


0.30 ±0.07 



tPNe coordinates from ICiardullo et al.l l|2004l ) . H II region IDs are from iBoulesteix et all 

lll974l 'l. with positions determined from ou r images. 

■t Galactocentric distance for D = 840 kpc jFreedman et ani200lh . 
*S = Subaru, K = Keck. 



Table 3. Reddening-corrected fluxes (blue). 



ID 


[On] 


[Nelll] 


H5 




H7 


[Oiii] 


He I 


Hen 


[Ariv] 


[Ariv] 




3727 


3868 


4101 


4340 


4363 


4471 


4686 


4711 


4740 


PN18 


36 ± 3 


91 ± 6 


26 


± 


2 


47 


± 3 


12.2 ± 0.8 


6.0 ± 3.4 


1.1 ± 0.3 


2.8 ± 0.3 


2.6 ± 0.3 


PN 55 




73 ± 5 


25 


± 


2 


47 


± 3 


8.5 ± 0.8 


5.1 ± 0.8 


(2.4) 


(2.4) 


(2.4) 


PN67 


145 ± 10 


24 ± 2 


25 


± 


2 


47 


± 3 


3.6 ± 0.5 


5.1 ± 0.7 


(0.5) 


(0.5) 


(0.5) 


PN68 


43 ± 8 


125 ± 9 


23 


± 


2 


47 


± 3 


26.9 ± 2.1 


(4.0) 


18.9 ± 1.4 


(2.6) 


5.2 ± 0.3 


PN 73 


71 ± 7 


66 ± 5 


27 


± 


2 


48 


± 3 


8.8 ± 0.9 


5.4 ± 0.6 


(1.5) 


(1.5) 


(15.5) 


PN 76t 


48 ± 4 


101 ± 6 


25 


± 


2 


47 


± 3 


14.3 ± 0.9 


5.0 ± 0.4 


1.3 ± 0.3 


2.1 ± 0.3 


2.5 ± 0.3 


PN 76* 


41 ± 3 


86 ± 6 


24 


± 


2 


47 


± 3 


13.3 ± 1.1 


4.6 ± 1.8 


2.0 ± 0.5 


2.8 ± 0.5 


3.0 ± 0.5 


PN 78t 


24 ± 14 


184 ± 13 


26 


± 


3 


47 


± 3 


58.7 ± 4.0 


8.8 ± 1.7 


58.7 ± 3.7 


13.3 ± 1.7 


14.8 ± 1.7 


PN 78* 


25 ± 4 


170 ± 12 


27 


± 


3 


47 


± 3 


59.3 ± 4.0 


7.3 ± 1.8 


62.9 ± 4.2 


13.3 ± 1.7 


15.9 ± 1.7 


PN83 


171 ± 15 


104 ± 7 


22 


± 


2 


47 


± 3 


20.7 ± 1.8 


(3.8) 


38.7 ± 2.4 


4.9 ± 0.8 


3.8 ± 0.7 


PN86 


47 ± 10 


113 ± 8 


26 


± 


2 


47 


± 3 


26.2 ± 2.0 


3.7 ± 1.2 


24.1 ± 1.6 


(2.2) 


(3.5) 


PN90 


69 ± 6 


68 ± 5 


23 


± 


2 


47 


± 3 


8.0 ± 0.8 


5.6 ± 0.9 


(2.7) 


(2.7) 


(2.7) 


PN94 


113 ± 9 


105 ± 7 


25 


± 


2 


47 


± 3 


22.8 ± 1.5 


(3.4) 


35.4 ± 2.3 


(3.2) 


(3.2) 


PN99 


35 ± 18 


130 ± 9 


27 


± 


3 


47 


± 3 


22.9 ± 2.1 


7.0 ± 1.6 


10.6 ± 1.3 


(3.2) 


6.8 ± 1.1 


PN105 


33 ± 5 


76 ± 5 


25 


± 


2 


47 


± 3 


10.9 ± 0.9 


(1.5) 


8.3 ± 0.7 


(1.5) 


(1.5) 


PN108 


117 ± 9 


130 ± 9 


26 


± 


2 


47 


± 3 


22.3 ± 1.9 


5.8 ± 1.2 


5.3 ± 1.2 


3.1 ± 1.1 


3.4 ± 1.1 


PNlll 






23 


± 


2 


47 


± 3 


20.0 ± 1.7 


6.4 ± 1.2 


9.7 ± 1.2 


6.2 ± 1.2 


6.4 ± 1.2 


PN113 


53 ± 4 


109 ± 8 


27 


± 


2 


47 


± 3 


15.5 ± 1.1 


6.0 ± 0.6 


8.3 ± 0.7 


2.2 ± 0.5 


2.6 ± 0.5 


B72 


187 ± 12 


5.8 ± 0.4 


23 


± 


1 


47 


± 3 


0.39 ± 0.04 


4.3 ± 0.3 


(0.1) 


0.3 ± 0.04 


(0.1) 


B90 


196 ± 12 


71 ± 4 


24 


± 


1 


47 


± 3 


5.2 ± 0.4 


6.7 ± 0.4 


4.4 ± 0.3 


(0.6) 


(0.6) 


B302 


155 ± 10 


6.1 ± 0.4 


21 


± 


1 


46 


± 3 


0.52 ± 0.04 


3.8 ± 0.2 


(0.1) 


0.3 ± 0.03 


(0.1) 



Values in brackets are upper limits. 
*Subaru/FOCAS spectrum. tKeck/LRIS spectrum. 
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Table 4. Reddening-corrected fluxes (red). 



ID 


[Oiii] 


He I 




Ha 




[Nil] 


He I 


[Sii] 


[Sii] 


[Ar III] 




5007 


5876 


6563 


6583 


6678 


6717 


6731 


7135 




1 ono _|_ vi 
izUo it /i 
























r^lN DO 


1 nc:o _|_ CO 
























pivr (i'T 


It Zo 
























rrs uo 


iZ40 it /O 


13.0 


± 


1.6 


276 


± 


32 


/I /I -1- 
44 it 


z.D it U.D 


(,2.0j 


z.y it u. i 


lZ. i it l.D 


P1\T 7*^ 


iZyO III iKj 


16.7 


± 


2.2 


303 


± 


39 


loo It iO 


o.o It u.y 




O.D It l.Z 


iU.O It z.o 


P1\T 7fit 
r 1 N i D ' 


iiOy it DO 


15.2 


± 


1.8 


282 


± 


33 


41 it 


o.o it U.o 


i.o It U.o 


Q -1- O /I 
O.O it U.4 


11 7-1-1 A 
ll. i it 1.4 




1 1 Q(i _|_ 70 
1±0D III /U 


























loyy it iuz 


22.4 


± 


2.8 


275 


± 


34 


c;7Q -1- 7n 
/ O it i U 


D.z it l.U 


10. i It z.u 


9K Q -I- q 
ZD.o it 0.4 


Zo.o it O.O 


PN78t 


1881 ± 113 
























PN 83 


1337 ± 79 


14.4 


± 


2.2 


281 


± 


33 


94 ± 11 


(3.1) 


22.5 it 2.9 


23.0 it 3.0 


19.6 it 2.5 


FN 86 


1734 ± 103 


13.2 


± 


1.9 


288 


± 


34 


53 ± 6 


(2.2) 


4.4 it 0.9 


7.3 it 1.1 


21.0 it 2.6 


PN 90 


943 ± 56 


14.6 


± 


1.8 


286 


± 


33 


43 ± 5 


4.2 ± 0.5 


11.8 it 1.4 


10.3 it 1.2 




PN 94 


1813 ± 107 


12.4 


± 


1.6 


280 


± 


33 


128 ± 15 


3.4 ± 0.7 


8.2 it 1.1 


9.6 it 1.2 




PN 99 


1593 ± 94 


21.2 


± 


2.8 


330 


± 


40 


45 ± 6 


5.8 ± 1.0 


5.5 it 1.0 


5.9 it 1.0 


24.1 it 3.1 


PN105 


1157 ± 68 


15.8 


± 


2.0 


279 


± 


34 


23 ± 3 


4.1 ± 0.6 


(1.2) 


2.5 it 0.5 




PN108 


1390 ± 83 
























PNlll 


1553 ± 92 
























PN113 


1418 ± 83 
























B72 


146 ± 9 


12.6 


± 


1.5 


309 


± 


38 


49 ± 6 


3.8 ± 0.5 


18.7 it 2.3 


13.7 it 1.7 


10.5 it 1.3 


B90 


698 ± 41 


14.3 


± 


1.7 


280 


± 


33 


42 ± 5 


3.5 ± 0.4 


34.2 it 4.0 


24.4 it 2.9 




B302 


182 ± 11 


12.6 


± 


1.5 


275 


± 


33 


33 ± 4 


3.5 ± 0.4 


13.4 it 1.6 


10.2 it 1.2 


9.9 it 1.2 



Values in brackets are upper limits. 
tSubaru/FOCAS spectrum. tKeck/LRIS spectrum. 
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Figure 5. Comparison between the emission line fluxes of vari- 
ous spectral lines, re lative to H/3, measured in this work and by 
iMagrini etlH | |2009|) . for 11 PNe in common. The line of one- 
to-one correspondence is flanked by lines drawn ±0.2 dex from 
it. Points representing O lines are shown in colour: [Oll]A3727 
(blue), [Olll]A5007 (green) and [Olll]A4363 (orange). The out- 
lier [O ll] A3727 line for PN 18 is identified. The red open squares 
represent the [Sll] AA6717, 6731 Unes. 



In principle, different excitation zones within the nebulae 
are characterized by different Te values, that would re- 
quire the detection of auroral lines of different ions, such 
as [Oll]A7325, [Slll]A6312 and [O ill] A4363. This is ob- 



servationally quite challenging for PNe at the distance of 
M33, and [O lll] A4363 is the only line that we are able to 
use to infer Te, with the exception of two [Nil] A5755 detec- 
tions (PN 76 and PN 78). Therefore, for the derivation of the 
ionic abundances we computed the line emissivities for all 
ions of interest at the electron temperature obtained from 
the [Oiii] A4363/(A4959 4-A5007) line ratio. This includes 
the low-excitation line [O ll] A3727, for which a different Te 
might be appropriate, but since most of the oxygen in bright 
PNe is doubly ionized, this introduces a negligible error in 
the final 0/H abundance. For PNe we did not use empiri- 
cal temperature schemes between high- and low-excitation 
zones, since, for reasons yet not fully understood, the un- 
certainties would be large (see, for example, the behaviour 
of Te derived from [N ll] A5755/(A6548 -t A6583) as a func- 
tion of Te from [ Qui] A4363/(A4959 -t A5007) for a sample of 
Galactic PNe in iGornv et al]|2009l ). For Hll regions, how- 
ever, we derived the temperature of the low-excitat ion zone 
for t he and N^ ionic abundances) from the iGarnettI 
1992h relation, T(0+) = 0.7 T(0++) + 3000 K, which is 
suppo rted by recent obse r vations of extragala ctic Hll re- 
gions (|Esteban et al.ll2009l : iBresolin et al.ll2009l ). 

The electron temperat ures were derived using the 
five-l evel program temden (|De Robertis. Dufour. &i HuntI 
Il987h. in the implem entation of iraf's nebular package 
( Shaw fc Dufouilll995l). The adopte d atomic parameters are 
the same as in Bresolin et al.l (|2009l ) , and are summarized in 
Table [5] for convenience. For Hel we used the emissivities of 
iPorter. Ferland. fc MacAdamI (|2007h . including the correc- 
tion for collisional excita tion, and for Hi and Hell those of 
iHummer fc StorevI (| 19871 ). 

The electron densities were estimated from the 
[S ll] A6717/A6731 ratios, when available (Subaru data only), 
or from [Ariv] A4711/A4740, if measured. We arbitrarily as- 
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Table 5. Sources of atomic data. 



Ion 


Transition 


Collision 




probabilities 


strengths 


[On] 


Frocse Fischer & 
Tachiev (2004) 


Taval (20071 


[Oiii] 


Froese Fischer & 
Tachiev (2004) 


Aegarwal & Keenan (1999) 


[Nil] 


Froese Fischer & 
Tachiev (2004) 


Hudson & Bell ('2005") 


[No III] 


Froese Fischer & 
Tachiev (2004) 


McLaughlin & Bell (20001 


[Ar III] 


Mendoza (1983) 
Kaufman Sz 
Sugar (1986) 


Galavis et al. (19951 


[Sll] 


Froese Fischer 
et al. (2006) 


Ramsbottom et al. (19961 


[Siii] 


Froese Fischer 
et al. (2006) 


Taval & GuDta (19991 



sumed A^e ~ 1000 cm ^"^ in the remaining cases. The density 
affects the Te determination only weakly, even though it has 
an important impact on the O"'' emissivity. For example, 
around 12,000 K, a variation of Ne between 1000 cm~'^ and 
5000 cm^'^ induces a change in T[Oni] of 60 K, which is on 
the order of 10% of the estimated temperature errors. The 
emissivity of the [O ll] A3727 line changes by about 50%, but 
we have verified that this has little effect on our results. For 
the final Te uncertainties (rounded to the nearest hundred 
K in Table [6l) we have propagated the 1-cr errors in the line 
fluxes (contained in Tables [3]|4]) , and include the small con- 
tribution from the uncertainties in A^^. Table [6] summarizes 
the plasma diagnostics thus obtained, and the density val- 
ues that we adopted. We note the generally good agreement 
between densities obtained from the [S ll] and [Ar iv] line ra- 
tios, when both are available. In these cases (PN76, PN 78 
and PN 83) we have adopted the [S n]-based densities for 
the following analysis, because of their smaller uncertainties. 
Table [6] reports the plasma diagnostics for the two targets 
in common between the Subaru and Keck datasets (PN 76 
and PN78). Within the errors, we find a good agreement 
between independent determinations for both the T[Oni] 
electron temperatures and the [Ariv]-based electron densi- 
ties. 

The abundances of 0+, 0++, Ne++, N+, Ar++ He+ 
and He"'"^, relative to H"^, were derived with the ionic pro- 
gram in IRAF, from the de-reddened line fluxes and the Te, 
Ne values found in the previous step. Total chemical abun- 
dances were o btained adopting the ionizati on correction fac- 
tors (ICFs) of lKingsburgh fc Barlow! (|l994h . with the excep- 
tion of argon for the Hll r e gions , in which case we adopted 
the recipes bv llzotov" et al.l(j2006l l. The ionization correction 
factor for oxygen is obtained as (1 -I- He^'''/He^)'^^^, where 
the He"*" abundance has been derived from He I A5876 when 
available (Subaru dataset), and from HelA4471 in the re- 
maining cases. Our final abundances are shown in Table [71 
where the uncertainties have been calculated by propagation 
of the errors in the line fluxes and Te- Once again, we point 
out the good agreement in the resulting abundances between 
independent measurements for PN 76 and PN 78, observed 



both with Subaru and Keck. The weighted abundances for 
these two PNe (also included in Table [7]) will be adopted for 
the remainder of the paper. 

For two of the Keck targets, PN 55 and PN 111, the 0+ 
abundance could not be calculated, because the [O ii] A3727 
line was outside the covered wavelength range, therefore the 
quoted value is a lower limit to the O/H ratio. However, 
the median 0^/0 ratio of the remaining sample is 0.06, 
exceeding 0.12 only in the case of PN 67 (0.26), so the total 
oxygen abundance in PN 55 and PN 111 is likely to be close 
to the values reported in Table [T] 

In Fig. [6] we compare our final PNe abundances with 
those obtained by M09. These authors have determined elec- 
tron temperatures either from direct detections of auroral 
hues ([Oin]A4363, [Nn]A5755), or from a statistical rela- 
tionship between Te and the strength of Hen A4686. Abun- 
dances obtained with the latter method are not expected to 
be accurate, neither on theoretical grounds, nor on the basis 
of the diagrams by M09 which were used to define this re- 
lation. In order to differentiate the two Te-determination 
methods, we represent the [O ill] A4363 detections in the 
M09 sample by full points in Fig. [51 while open symbols are 
used for PNe without a A4363 line detection (we remind the 
reader that [O lll] A4363 is measured for each PN contained 
in our sample). It can be seen that the agreement between 
our abundances and those by M09 is generally acceptable, 
and that no obvious systematic offset exists between the 
two sets of abundances, with the exception of helium. This 
is likely a result of different abundance calculations, in par- 
ticular regarding the effect of coUisional excitation, because 
the He I line fluxes are in good agreement between the two 
studies. The remaining large discrepancies (> 0.2 dex) in 
Fig. [6l are due to significant differences either in the mea- 
sured diagnostic line strengths (as in the case of Ne), or in 
the adopted Te (in the case of the most discrepant points for 
O, which do not have a [O in] A4633-based Te determination 
in M09). The plot for O confirms the expectation that PNe 
in the M09 sample without an [O ill] A4363 detection, and 
whose abundances derive from a Te determination based on 
He II A4686, are affected by larger random errors than those 
with an observed [O ill] A4363 line, and that a large abun- 
dance scatter is expected for these objects. 

One of the three Hll regions we observed, 
B 90 , was also in cluded in the abundance studies 
by ICrockett etlo] (|2006| 1 and RS08. Both obtained 
12 -I- log(0/H) = 8.50±0.06, in very good agreement with 
our value 12 -|- log(0/H) = 8.47±0.08. The Ne/H abundance 
ratios are also virtually coincident. 

For the remainder of the paper, we complement our 
dataset of 16 PNe, covering mostly the inner region of M33, 
with the PNe sample studied by M09, but considering only 
objects with a direct determination of the oxygen abun- 
dance, i.e. based on the measurement of the [O ill] A4363 
auroral line. This restricts the number of PNe from M09 to 
32 (out of 93), extending in galactocentric distances from 
0.6 to 8 kpc. Of these, 6 are in common with our sample: 
PN18, PN67, PN73, PN76, PN108, and PN113. New au- 
roral line-abundances for 10 PNe are therefore added by our 
study, 7 of them at 7? < 1.2 kpc. For consistency, we have 
recalculated the chemical abundances from the line fluxes 
published by M09, with the same procedure used for our 
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Table 6. Plasma diagnostics. 



ID 


T[0 III] 


T[Nii] 


JV,([Sii]) 


Afe{[Ariv]) 


(adopted) 


PN 18 








Qoqn _|_ eon 


Qonn -1- 890 


PN 55 


1 0700 4- '^OO 








1 noo -1- 1 00 


PN 67 


1 0400 4- fiOO 








1 noo 4- 1 00 


PN 68 


1 '^TOO 4- 1 000 








1 noo 4- 1 00 


PN 73 


1 ni 00 4- '^OO 




1 fif^O -1- 1 90 




1 n«n 4- 1 90 
J-UUU ^ ±zu 


PN 76^ 


1 9'^00 4- fiOO 




finoo -1- 1 490 


QQon -1- 1 ^^90 


«n90 4- 1 490 


PN 76^ 


1 1 nnn -U fiOO 










PN 78"'' 


ofil 00 -k 1 fiOO 




4Q1 n _|_ Qon 






PN 78* 


18900 ± 1400 






9360 lb 1310 


9360 ± 1310 


PN83 


13500 ± 800 




540 ± 50 


1100 lb 990 


540 ± 50 


PN86 


13400 ± 800 




2980 ± 840 




2980 ± 840 


PN90 


10900 ± 600 




240 ± 20 




240 ± 20 


PN94 


12500 ± 600 




850 ± 80 




850 ± 80 


PN99 


13200 ± 800 




660 ± 150 




660 ± 150 


PN105 


11200 ± 500 








1000 ± 100 


PN108 


13600 ± 900 






6310 lb 3360 


6310 ± 3360 


PNlll 


12500 ± 700 






5230 lb 1530 


5230 ± 1530 


PN113 


11700 ± 600 






7180 lb 2150 


7180 ± 2150 


B 72 


7800 ± 300 




51 ± 5 




51 ± 5 


B 90 


10400 ± 400 


10200 ± 1400 


22 ± 2 




22 ± 2 


B 302 


8000 ± 300 


8100 ± 700 


79 ± 8 




79 ± 8 



tSubaru/FOCAS spectrum. tKeck/LRIS spectrum. 



PNe sample. In the plots that follow, in the cases where 
the M09 and our sample overlap, we have retained only the 
abundance information obtained from our observations, in 
order to avoid representing the same object more than once. 



As mentioned above, ICiardullo et al.l (|2004l ) identified 
two PNe (PN24 and PN67) as possibly belonging to the 
M33 halo, having velocities that deviate more than 2.5 a 
from the mean rotation curve. However, they point out 
that these outliers can be expected out of the total num- 
ber of 140 targets, simply on statistical grounds. We find 
that the chemical abundances of PN 67 are comparable with 
those of the other targets in our sample. PNe in the halo 
of the Milky Way ar e generally metal-poor compared to 



those in the disk (e.g. Howard. Henrv. fc McCartney! [l997l : 



iHenrv. Kwitter. fc Bailed 120041 '). We conclude that PN67 



likely belongs to the disk of M33, rather than the halo. 
However, we also point out that PN 67 is the only object in 
our PN sample that clearly shows absorption components 
in the Balmer lines, an indication of the presence of an un- 
derlying stellar component, as in the case of H ll regions. It 
has also the lowest excitation ratio observed in our sample, 
O++/O = 0.74, which overlaps with the range observed for 
H 11 regions in M33 (however, PNe wit h even lower excitation 
are k no wn, e.g. in the Mil ky Way, see lKingsburgh fc Barlowl 
Il994l or lGornv et~alll2009). PN 67 is not among the brightest 
in Hq flux or in msoor ( CiarduUo et al.ll2004h . It is there- 
fore unclear whether this is truly a PN, or a compact, high- 
excitation H II region. For the remainder of this paper we will 
continue to include PN 67 along the remaining PNe. None 
of the results presented in the next sections suggests that 
we should do otherwise. 



cn 
o 



cn 

o 

+ 

CN 




7.6 8 6 6.4 

12 + log(X/H) (this work) 

Figure 6. Comparison between the chemic al abundances of He , 
O, Ne and Ar measured in this work and bv lMaerini et aL 1 1|2009| . 
=M09) for PNe in common. Dashed lines are drawn to represent 
the one-to-one correspondence. Objects with an electron temper- 
ature determination based on [O III] A4363 from the M09 sample 
are shown as red full dots. 



4 EXCITATION PROPERTIES 

A flrst look at the main properties of the PNe in M33 
is presented in Fig. [T] where we show the He^"'"/H and 
N/O abundance ratios as a function of the excitation 
0^^/(0^ -I- O^^) (panels a and b). In this and subsequent 
plots we use red open circles to represent our new PN sam- 
ple, while the green circles represent the PNe observed by 
M09 and with a [O ill] A4363 line detection. We also show 



PNe m M33 11 



Table 7. Chemical abundances. 



ID He/H 12 + log(0/H) 12 + log(N/H) 

PN18 0.12 ± 0.07 8.47 ± 0.09 

PN55 0.10 ± 0.02 8.49*± 0.10 

PN67 0.10 ± 0.02 8.31 ± 0.12 

FN 68 0.11 ± 0.01 8.16 ± 0.10 8.07 ± 0.11 

PN73 0.12 ± 0.02 8.69 ± 0.10 8.62 ± 0.12 

PN76t 0.10 ± 0.01 8.38 ± 0.09 7.89 ± 0.10 

FN 76* 0.09 ± 0.04 8.43 ± 0.10 

PN76+ 0.10 ± 0.01 8.40 ± 0.09 7.89 ± 0.10 

PN78t 0.17 ± 0.02 8.11 ± 0.10 9.32 ± 0.11 

FN 78* 0.18 ± 0.03 8.21 ± 0.10 

PN78+ 0.17 ± 0.01 8.15 ± 0.07 9.32 ± 0.11 

PN83 0.14 ± 0.02 8.41 ± 0.10 8.01 ± 0.11 

FN 86 0.11 ± 0.01 8.49 ± 0.10 8.27 ± 0.11 

PN90 0.11 ± 0.02 8.45 ± 0.10 7.99 ± 0.11 

PN94 0.12 ± 0.01 8.63 ± 0.09 8.47 ± 0.10 

FN 99 0.17 ± 0.02 8.42 ± 0.10 8.36 ± 0.11 

FN 105 0.12 ± 0.02 8.49 ± 0.09 8.04 ± 0.11 

FN 108 0.11 ± 0.02 8.36 ± 0.11 

FN 111 0.13 ± 0.02 8.47*± 0.10 

FN 113 0.12 ± 0.01 8.54 ± 0.09 

B72 0.09 ± 0.01 8.46 ± 0.09 7.45 ± 0.09 

B90 0.11 ± 0.01 8.47 ± 0.08 7.53 ± 0.08 

B302 0.09 ± 0.01 8.44 ± 0.08 7.36 ± 0.08 

* Lower limit. 

tSubaru/FOCAS spectrum. tKeck/LRIS spectrum. 

+Adopted abundances (weighted mean of Subaru and Keck results) 



12 + log(Ne/H) 12 + log(Ar/H) 



7.79 ± 0.10 
7.79 ± 0.11 
7.48 ± 0.13 
7.54 ± 0.11 
7.87 ± 0.12 
7.75 ± 0.10 
7.72 ± 0.11 
7.74 ± 0.10 



7.48 
7.51 
7.50 
7.71 
7.71 
7.77 
7.81 
7.74 
7.76 
7.74 



± 
± 
± 
± 
± 
± 
± 
± 
± 
± 

7.74 ± 0.12 
7.86 ± 0.11 
7.62 ± 0.12 
7.94 ± 0.10 
7.52 ± 0.10 



0.11 
0.11 
0.09 
0.11 
0.11 
0.12 
0.10 
0.12 
0.11 
0.12 



95 


± 





10 


44 


± 





11 


11 


± 





09 


11 


± 





09 


13 


± 





11 


13 


± 





11 


25 


± 





11 


29 


± 





10 



6.36 ± 0.11 

6.31 ± 0.11 
6.25 ± 0.10 



with orange squares our new Hll region abundances (B 72, 
B 90 and B 302), and with the remaining squares a sample of 
Hll regions for which auroral line-based abundances could 
be derived, using line fluxes we compiled from the literature 
("Kwittc r fc AUerl Il98ll: IVflchez et all Il988l : ICrockett et"al] 
[2OO6: Magrini et al.ll2007 : RS08). For consistency, we recom- 
puted the PN and H 11 region abundances adopting the same 
techniques explained in Sect. 3. 



Fig. [7] highlights the fact that the excitation of PNe is 
very high, with more than 75% of the oxygen in the O^^ ion- 
ization stage for most of the objects shown. The Hll regions 
can be clearly separated from the majority of the PNe in 
Fig- EK simply from their He^"'"/H ratios, since HellA4686, 
when observed in H 11 regions, never exceeds the intensity of 
a few percent of H/?, while in planetary nebulae this line is 
often observed and can reach intensities similar to H/3. The 
separation between PNe and Hll regions is well displayed 
also by looking at N/O vs. 0++/(0+ + 0++) in Fig.[7j3. The 
average N/O ratio of the PNe is obviously larger than that 
for the H 11 regions, a result of the nucleosynthesis processes 
taking place in the advanced stages of evolution of the PNe 
stellar progenitors. This is also reflected by the large dis- 
persion seen for the PN data, compared to the Hll regions. 
Finally, in Fig. [T]; we plot the ionic temperature Te [O ill] , 
derived from the [Olll] A4363/(A4959 -f A5007) line ratio, as 
a function of the oxygen abundance. The general trend, for 
both PNe and Hll regions, is of decreasing electron tem- 
perature with increasing metallicity. However, at any given 
0/H ratio, the PNe have higher Te values, as expected, and 
display a larger dispersion, which reflects the dispersion in 
stellar effective temperatures. 



-0.5 



O 



O 




0.4 0.6 
0+V(0^ + 0++) 

Figure 8. The Ne/O ratio as a function of the excitation ratio 
0++/(0+-|-0++). Symbols as in Fig. [7] The error bars for the 
H II regions are shown. The hori zontal dashed line re presents the 
solar value log(Nc/0)o = -0.76 llAsplund et al.ll2009l ). The mean 
Ne/O = 0.20 for the PNe in our sample is shown by the continuous 
blue line. 



Uncertainties in the ICFs? 

Fig. [8] shows how the Hll region Ne/O ratio appears to 
decrease with excitation, in contrast with the expectation 
that this ratio should remain constant with variations of the 
nebular ionization structure. Since our neon abundances are 
derived from the single [Nelll]A3868 line, this result can 
originate from an inaccurate correction for unseen ioniza- 
tion stages, in particular Ne"^, when the 0^'*"/(0^-|- O^^) 
ratio is below 0.5. However, it is unclear whether the dis- 
crepancy at 0'^'''/(0"'""^-|-0'^) < 0.5 is due to problems with 
the ICF(Ne) or to systematic errors in the data, because 



12 F. Bresolin et al. 



PNe 

O this paper 

o Magrini et al. (2009) 



X 0.05 



0) 

X 







0^ 



0.5 
V(0+ + 0++ 



) 



H II regions 

■ this paper 

□ Crockett et al. (2006) 

□ Kwitter & Alter (1981) 



1 1 1 1 1 1 1 1 
O o 




o ° 


O 
\ 


o o - 




- 






o 






— on] 1 ■■■Hill 1 [nnnanKJ®o — 

1 1 1 1 1 1 1 1 





-III 


1 1 1 1 j_ 




b ~ 




oo o - 






-nu 

1 1 1 1 1 1 




— 1 1 1 


1 , , , 1- 



0^ 



0.5 

V(o+ + 0++ 



) 



□ Vilchez et al. (1988) 

□ Magrini et al. (2007) 

o Rosolowsky & Simon (2008) 



20 - 



O 

o" 10 




8.5 9 
2 + log(0/H) 
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legend at the top. 



some of the low-excitation objects have a 'normal', approx- 
imately so lar-like Ne/0 ratio [log (Ne/O)© = -0.76 ± 0.11, 
lAsplund~et al.., 20091 . shown by the horizontal dashed line]. 
Moreover, some of the high-excitation Hll regions signifi- 
cantly exceed both the solar value and the PN values. We 
have considered alternative scheme s for the ICF ( Ne) o f H 11 
regions, such as thos e prop osed by llzotov et al.1 (|2006l ) and 
IPerez-Montero et al.l (|2007l ). but found even larger trends 
with excitatio n. We have therefore c o ntinu ed to adopt the 
ICF(Ne) from iKingsburgh fc Barlow! (|l994 ) for the rest of 
this paper, in order to minimize the excitation dependence 
of the Ne/0 ratio. 

Considering only our PN sample, we obtain a mean 
<Ne/O> = 0.20 ± 0.03, shown by the continuous line in 
Fig. [HI We suggest to use this value as the representative 
Ne/0 ratio in the ISM of M33, because the selected PNe 
are the objects with the highest excitation (and presumably 
with the most accurate abundances), and therefore the un- 
certainty in the total Ne abundance due to corrections for 
Ne"*" should be negligible. 

As we explained earlier, for a r gon w e have adopted the 
ICFs from |Kingsburgh fc Barlow! (| 19941 ) for PNe and from 
!lzotov et "alT i 2006! ) for H 11 regions. A plot of the Ar/O abun- 
dance as a function of the excitation 0'''^/(0^''' -I- O"*") is 
shown in Fig. 19] As in the case of neon, the Hll regions show 
a trend, this time with high-excitation objects having a lower 
Ar/O ratio than lower- excitation ones. The la tter agree with 
the solar ratio, log(Ar/O)0 = -2.14 ± 0.09 (!Loddersl!2o"o^ . 
indicated by the horizontal dashed line in Fig.[9|. This sug- 
gests that the ICF(Ar) we are adopting is inadequate at 
0++/(0++ -I- 0+) > 0.6. We also point out the existence of 
a systematic offset of the Ar/O ratio between our PN sample 
(red circles), which agrees with the solar value, and the one 
based on the M09 data, which has a lower mean Ar/O. This 
could be related to observational errors, or could be due to 
the systematically higher extinction found by M09 (which 
would lower the [Arm] A7135 strength relative to H/3). For 
the following discussion we must therefore keep in mind that 
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high-excitation H 11 regions and the M09 sample of PNe ap- 
pear to have argon abundances that are under-estimated. 



5 ABUNDANCE PATTERNS 
5.1 Nitrogen and helium 

The impact of the nucleosynthetic processes occurring in 
the PN progenitors on the chemical compositions of the 
ejected envelopes can be visualized in a plot of N/O vs. 
He/H (Fig. nU)) . Both N and He can be enhanced at the 
stellar surface during the red giant and AGB phases, largely 
depending on the mass and metallicity of the stellar precur- 
sors. The details are complicated by the combined effects of 
mechanisms such as ma ss loss, rotation, h ot bottom burn- 
ing and third dredge-up (!Charbonnel!!2005l ). The comparison 
between our new PN sample (mean log(N/0) = — 0.27, ex- 
cluding PN 78) and the Hll region data (square symbols, 
with a mean log(N/0) = —1.22) shows a N/O enhancement 



PNe m MS 3 13 



-2 - 



PN 78 
) 1 



PN 99 




_I I I L_ 



_I I I L_ 



0.05 



0.15 



0.2 



He/H 



Figure 10. The N/O abundance r atio as a function of the He/H 
ratio. The boundaries defined by iTorres-Peimbert fc Peimber j 
1I1997I ') for Type I PNe in the LMC (N/O > 0.38, He/H > 0.105) 
are drawn in orange. Symbols as in Fig. [7] Error bars are shown 
for our sample. The mean N/O ratio is drawn with dashed hor- 
izontal lines for the H II regions (lower line) and our PN sample 
(upper line). 



on the order of one dex (the two mean values are shown as 
dashed lines in Fig. I10|l . When we factor in the radial trends 
(Sec. I6.2]l . the typical enhancement in N/O is ~ 0.8 dex. 

Stars belonging to the high end of the PN progenitor 
mass spectrum (which extends approxima tely between 0.8 
and 8 Mq) are associated with Type I PNe (jPeimbc rt,, 197^ . 
and are part of the youngest stellar population, showing 
clear N and He abundance excesses relative to the compo- 
sition of H II regions, which represent the chemical makeup 
of the present-day ISM. The values of the minimum N/O 
and He/H ratios that define this class of PNe vary slightly 
between different authors, and need to account for the value 
of the N and O abundances in the ISM out of which the pro- 
genitor stars were formed (e.g. IStasihska. Richer, fc McCalll 



Il998l : iLeisv fc Dennefeldlll996l l. Since the metal content of 
M33 is comparable to that of the LMC, we show in Fig. (TO] 
the corr esponding limits, N/O > 0.38 , He/H > 0.105, pro- 
posed bv lTorres-Peimbert fc Peimber^ (jl997^ . With this def- 
inition, most of the PNe in our sample (red symbols with 
error bars) are of Type I, or have an N/O ratio that is 
compatible, given the errors, with being of Type I. How- 
ever, it is worth pointing out that most of the PNe in our 
sample have N/O ratios that are only slightly above the 
iTorres-Peimbert fc PeimbertI l)l997h limit for inclusion into 
the Type I definition. In fact, we do not expect that a large 
number of PN progenitors in our sample, which is drawn 
from the brightest PNe in M33, experienced evolutionary 
phases that can strongly enhance the N/O ratio, such as 
hot bottom burning or second dredge-up, because the cor- 
responding progenitor masses would b e arou nd 2.5 Mq, ac- 
cording to the models of iMarigo et al.l (|2004l ) . Furthermore, 
the IMF would favor the detection of lower mass progenitors 
with respect to high mass ones. 



5.1.1 The remarkable PN 78 

The largest helium content in our PN sample (He/H ~ 0.17) 
is reached for PN 78 and PN 99 (both are identified in 
Fig. I10|l . The former also has a striking nitrogen excess. 



12 + log(N/H) = 9.20 and log(N/0) = 1.09. The large N con- 
tent does not depend on the correction factor scheme that 
we used, and is due to the strong nitrogen lines observed 
(the strength of [N 11] A6583 is approximately twice that of 
Hq; see the spectrum in Fig. [2]). The high electron tempera- 
ture is accompanied by a very high nebular excitation, with 
the only detection of [Ar v] AA6435, 7005 in our sample. 

The nitrogen enrichment that we derive for PN 78 is 
comparable to the maximum obser ved to date in Type I PN e 
belonging to the Milky Way (e.g. IStanghellini et al.1 [iooil ) . 
although it occurs at a much lower O/H ratio. In fact, PN 78 
has the lowest measured oxygen abundance in our sample, 
12 -I- log(0/H) = 8.15, considerably lower than the median of 
8.47. It is tempting to interpret these data as the result of 
hot bottom burning in a massive progenitor (~ 5 Mq), which 
could explain both the low O and the high N abundances, 
with nitrogen synthesized at the expense of C. At sufficiently 
high temperatures N could be produced at the expense of 
oxygen in the ON cycle. However, in the case of PN 78 the 
latter possibility seems to be excluded. We plot in Fig. llll the 
N/O ratio as a function of N/H. The regression line to the 
whole PN sample (excluding PN78), calculated accounting 
for errors in both co ordinates (we use d an IDL adaptation of 
the routine fitexy bv lPress et al.ll99^ for this and subsequent 
fits of this kind), has a slope of 0.99±0.01, indicating that the 
N/O enhancement relative to H 11 regions is due primarily to 
high N production, rather than O depletion. This appears 
to hold for the M33 PN sample as a whole, and also for the 
specific case of PN 78. In addition, both the Ar/O and Ne/O 
ratios for PN 78 are comparable to those of the rest of the 
sample, confirming that oxygen is not significantly depleted. 

PNe observed in the Milky Way and the Magellanic 
Clouds have generally much smaller N/H ratios than we 
find in PN78, which is interpreted as a modest efficiency 
of hot bottom bu rning in intermediate-mass precursors 
(jMarigo et al.ll2003l ). since a larger efficiency would lead to 
an overproduction of N. Our data for PN 78 seem to suggest 
that, occasionally, large N/H and N/O ratios can be encoun- 
tered, simultaneously with a large He/H ratio (the latter as 
a re sult of an efficient third dredge- up). 

[Richer fc McCalll (|2007l ) have presented chemical abun- 
dances for PNe in the irregular galaxy NGC 6822. 
One of their targets, S33, with 12 + log(N/H) = 8.93 and 
log(N/0) = 0.90, has a nitrogen excess similar to PN78 
in M33. In this case, however, the oxygen abundance, 
12 -I- log(0/H) = 8.03, is Zan;er than the galaxy's mean value. 
These authors propose that this could be an example of a 
star in which oxygen has been dredged-up in the late evo- 
lutionary plujseS;_Ilie_Jugh N/O ratio is not explicitly dis- 
cussed bv lRicher fc McCalll (|2007f l. although they note that, 
in general, nitrogen should be more easily dredged-up than 
oxygen. 



5.2 Neon and oxygen 

Since core-collapse supernovae (massive stars with initial 
AI > 8Mq) are the main sites of oxygen and neon pro- 
duction, one expects that, unless important nucleosynthesis 
of these elements takes place in the less massive stars that 
are the PN precursors, their variation occurs in lockstep in 
PNe, similarly to the case of Hll regions. This is in f act ob- 
served in the Milky Way and the Magellanic Clouds jHenrvl 
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Figure 11. The N/O abundance ratio as a function of tlic N/H 
ratio. Symbols as in Fig. [T] The regression to the PN data is 
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1 19891 : iLeisv fc Dennefeldl |2006| ). and is normally taken as 
evidence that the oxygen abundance measured in PNe is 
equivalen t to that of the progen itor stars at the moment 
of birth (|Richer fc McCalll 120071 ). Deviations from this re- 
sult have been discovered in low metallicity environments, 
where third dredge-up episodes could be more efficient, with 
the consequence that oxygen can be overabundant in PNe 
relative to H 11 regions, as b oth th eorv (|Marigdl2001^ and ob- 



servations dMagrini et aL 2005'; Pc na. Stasiiiska. fc Richer] 
l2007l :lKr 



according to [Richer 



Cniazev. Pustilnik. fc Zuckcr, 20081 ) suggest. However. 

McCalll (|2007l ) these findings consti- 
tute the exception rather than the rule, perhaps because 
most progenitors of the bright PNe studied spectroscopi- 
cally in other galaxies are of rather small mass (< 1.5 M©, 
[Richer fc McCallll2008l ). 

Fig. 1 121 shows the values of Ne/H as a function of 0/H 
for PNe and Hll regions, with the same symbols as in Fig.[71 
The regression to the PN data (considering [O ill] A4363- 
based abundances only), shown by the dashed line, has a 
slope of 1.09 ± 0.07, with a correl ation coefficient of 0.89. 
The Ne-0 relationship observed bv llzotov et al.l (|2006l ) in a 
sample of star-forming galaxies and blue compact galaxies 
is shown for reference (continuous line, slope 1.097 ± 0.015). 
The increasing deviation of the RS08 H 11 region data (blue 
squares) from the Izotov et al. line with decreasing O/H has 
the same nature as the excitation trend we detected in the 
Ne/O ratio (Fig. [Hjl, and can be related to an inaccurate 
ICF (Ne) at low n e bular excitation. 

IWang fc Liul (|2008l ) have pointed out that the Ne/O 
ratio of the ionized ISM follows an increasing trend with 
O/H, reaching a mean value of ~0.25 in the Milky Way and 
other metal-rich galaxies, suggesting that the Ne enrichment 
of the ISM is delayed relative to O. The average Ne/O ratio 
for PNe with [O ill] A4363 detections in M33 (our sample -f 
M09) is 0. 18 ±0.04, com parable with the value observed in 
the LMC (|Leisv fc Denn cfcld 200i) , which has a present-day 
metallicity similar to that of the ISM in M33. 



5.3 Argon and oxygen 

Argon and sulfur are also products of primary nucleosyn- 
thesis in massive stars, and their abundances can be used in 
alternative to oxygen to characterize the initial composition 
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Figure 12. No/H as a function of O/H. Symbols as in Fig. [7] 
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Figure 13. The Ar/H abundance ratio as a function of O/H. 
Symbols as in Fig. [7] The regression to the PNe (our sample + 
M09) is shown by the short-dashed line, and the trend found by 
I Izotov et al. ] 1 I2OO6I) for more than 400 emission-line galaxies is 
shown by the continuous line. The long-dashed line represents 
the regression to our PN sample, combined with the H II regions. 



of PN progenitors, since they s hould not be affected by evo- 
lution during the AGB phase (|Leisv fc Dennefeldll2006h . at 
any chemical composition. The main difficulty in the study 
of argon and sulfur is that observationally it is more chal- 
lenging than that of oxygen, due either to the weakness of 
the lines, or the required wavelength coverage. Moreover, 
the correction for unseen s tages of ionization can be impor- 
tant, especially for sulfur (|Henrv et al.ll2003 ). In this work 
we only consider the abundance of argon, since for sulfur we 
only have access to the [S 11] lines. 

Observations in the Milky Way have shown that argon, 
like neon, fo llows oxygen reaso nably well in both PNe and 
H II re gions (|Henrv et al.ll2003 ). although IStanghellini et"al] 
(|2006h failed to find a correlation between Ar and O abun- 
dances for PNe, likely as a result of the uncertain Ar abun- 
dances. Our plot of Ar/H vs. O/H is shown in Fig. 1131 where 
we can see a good correlation {Rxy ~ 0.76) between the two 
quantities. The slope of the regression that considers the 
whole PN sample (short-dashed line) is 0.58 ± 0.09, signif- 
icantly lower than unity, and smaller than the value found 
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in kpc (top axis). Symbols as in Fig.[7l The dashed line represents 
the linear regression to the H II region data. 



for th e star-forming galaxy sample studied by 
20061) ■ 0.916 ± 0.021. For Milky Way PNe, 



Izotov et al 



Henry et al 



20041 ) obtained a good correlation between Ar/H and O/H, 



but with an even higher slope of 1.34 ± 0.14. In light of the 
problems with the PN Ar/H abundances mentioned earlier, 
we carried out a regression to our PNe sample and the H 11 
regions combined, since they appear to follow a common re- 
lation in the diagram (long-dashed line). We find a slope of 
0.98 ±0.27, consistent with unity, although the error is quite 
lar ge. We point out th at these data lie systematically above 
the llzotov et all (|2006l ) regression line. We conclude that Ar 



is rather well correlated with O in both PNe and H 11 regions 
in M33, but further studies will be required to reduce the 
uncertainty in the relation between these two quantities for 
this galaxy. 



6 RADIAL ABUNDANCE GRADIENTS 

One of our main motivations at the inception of this spectro- 
scopic study of M33 was to test whether PNe could provide 
complementary information on the galactic-scale oxygen 
abundance gradient in this galaxy and on its possible tempo- 
ral evolution. Despite the fact that several Hll region emis- 
sion line studies have been carried out in M33 throughout 
the ye ars (to cite a few:lKwitter fc Alleilll98ll:IVflchez et al.l 



198^: IWiUner fc Nelson-Pate]H2002l: ICrockett et al.1 I2OO6I: 



Magrini et al.ll2007l : lRosolowskv fc Simonll2008l : lRubin et all 
20081 ). the present day metal content of the central parts of 



the galaxy remains the most uncertain, due to the difficulties 
in measuring the required, temperature-sensitive faint auro- 
ral lines, as the gas temperature decreases with decreasing 
galactocentric distance. This is a consequence of the neb- 
ular cooling increasing with metallicity. Fig. [14] illustrates 
the systematic change of temperature with distance from 
the galactic center derived for the Hll regions, for which 
we find a gradient of 300 ± 60 Kkpc~^ (no galactocentric 
variation in temperature is detected for PNe). 

Comparative studies of the abundance gradient in the 
Milky Way, that include observations of PNe, Hll regions 
and young stars, provide rather consistent results among 
the different abu ndance tracers (jMaciel fc Quirezal 1 19991 : 
iHenrv et al.ll2004l ). but with considerable uncertainties re- 



sulting from poorly known distances to individual targets. 
Chemical elements that are unaffected by nucleosynthesis 
in low- and intermediate-mass stars, such as S, Ar and CI, 
should preferentially be used, however it appears that O and 
Ne, which are more easily measured, also trace the abun- 
dances of the ISM at the time of formation of the PNe pro- 
genitors, since their production is very small in the corre- 
sponding mass range, except at the lowest metallicities. 

Our discussion about the Ne/O ratio presented in Sec- 
tion [52] confirms the general finding obtained by other au- 
thors in a few nearby galaxies, namely that above the metal- 
licity of the SMC [12-Hog(0/H) =8.1] the brightest PNe 
and H 11 regions h ave comparable metallicities, as me asured 
by the O/H ratio (|Richeij|l99 AlRicher fc McCallll2007l ). This 
allows us to examine the metallicity gradient in M33 inde- 
pendently of the Hll regions, by using PNe. In the compar- 
ison, we can also include massive supergiant stars, which, 
owing to their young ages, are also excellent probes of the 
present day chemical composition of galaxies. In a recent 
work on the spiral galaxy NGC 300, iBresolin et~aD (j2009l ) 
confirmed that the analysis of the metal lines of blue super- 
giants provides an abundance gradient which is in excellent 
agreement with the one derived from [O in] A4363 detections 
in Hll regions, including metallicities near the solar value. 
In the case of M3 3, stellar metalli c ities f or young stars have 
been published bv lUrbaneia et all (|2005l . B supergiants) and 
lU et al.1 (|2009l . A supergiants). The latter also compared the 
radial abundance gradient in M33 using different indicators 
(Hll regions, supergiant stars and Cepheids). 

Despite the difficulties of obtaining spectra of extra- 
galactic PNe with good signal-to-noise ratios, the use of PNe 
presents a couple of advantages over Hll regions in abun- 
dance studies near the galactic central regions. First, as the 
metallicity increases at small galactocentric distances, the 
determination of the decreasing Hll region electron tem- 
perature becomes very difficult, as a consequence of the 
more efficient gas cooling that reduces the strengths of the 
[O III] lines. In the case of PNe, the excitation provided 
by the much harder stellar radiation and the higher elec- 
tron densities, which reduce the cooling due to coUisional 
de-excitation, ensure that the detection of [O ill] A4363 is 
still feasible. Secondly, PNe are not expec ted to be affected 
by abundance biases at high O/H ratios. IStasihskal (|2005l ) 
pointed out that in the presence of temperature gradients 
within Hll regions, which are predicted to develop at high 
metallicity (around and above the solar value), the mea- 
sured O/H abundances could systematically underestimate 
the true metal content of the gas. Again, the high effective 
temperature of the ionizing stars of PNe, together with their 
higher densities, prevents this from occurring. 

In the following discussion, we consider optical 
comparison samples of H 1 1 regio ns taken from vari- 



ous authors: iKwitter fc AUeil Il98l]). IViIchez et all (|l988l ). 



[Crockett etHl (|2006l ). [Magrini et all (|2007[ ) and RS08 



These works provide [O ill] A4363 detections for the measure- 
ment of the electron temperature. The wavelength cover- 
age differs among the different studies, so th at, for example, 
[Nil] a nd [Ar ill] fluxes are unavailable for the [ Crockett et al] 
(|2006l ) and RS08 samples. Only nebulae for which the elec- 
tron temperature could be directly derived from [O ill] A4363 
were retained, leaving a flnal sample of 88 nebulae. To 



16 F. Bresolin et al. 



R (kpc) 



8 



9 - 



> 8.5 = OQ 



X 

\ 

cn 
O 

+ 

CNI 



8 - 



7.5 



94 O 



□ o 

O 7J 



O q--— ------ "--io. 



□ O □ 



□ ^ 

o o 

^ 

o □ 



SB OO 78 



PNe 

O this paper 

O Magrlni et al. (2009) 

_^ ^ ^ \ L_ 



H II regions 

□ this paper □ Vilctiez et al. (1988) 

□ Crocl<ett et al. (2006) □ Magrini et al. (2007) 

□ Kwitter & Alter (1981) « Rosolowsky & Simon (2008). 

J ^ I ^ \ ^ ^ ^ ^ ^ ^ ^ 







0.2 



0.4 0.6 
R/R25 



0.8 



1 



Figure 15. O/H radial gradient for PNe and Hll regions. Radial distances are given both in terms of the isophotal radius R25 (bottom 
axis) and in kpc (top axis). The different symbols used are summarized in the legend. Linear regressions are shown by the red continuous 
line (PNe) and the dashed black line (Hll regions). The four planetaries that deviate the most from the PNe regression, PN68, 73, 78 
and 94, are labeled. 



overcome difficulties in the analysis arising from the non- 
homogeneous nature of the data sets and to avoid systematic 
discrepancies that could originate from differences in atomic 
data, ionization correction factors, etc., we have recalculated 
the abundances from the published line fluxes using the same 
procedure employed for our Hll region data (see Section |3]). 

6.1 Oxygen 

The O/H radial behaviour for the PNe having [O ill] A4363- 
based abundances is compared with the H 11 region and the 
blue supergiant gradients in Fig. ll5l and ll6l respectively. We 
have calculated a linear regression for the combined PNe 
data (ours and M09), and represent it with the continuous 
red lines in both figures. W e have used the weighted l east 
square method described by lAkritas fc Bershadvl (| 19961 ) , to 
account for the intrinsic scatter in abundance, in addition 
to uncertainties in the individual abundance values. The 
regression yields: 

12 + log(0/H) = 8.47 (±0.07) - 0.013 (±0.016) Rkpc (1) 

The gradient slope that we obtain is compatible within 1-a 
with both a flat gradient and with the slope obtained by 
M09 (-0.031 ±0.013 dexkpc"^). We note that M09 did not 
restrict their analysis to objects with auroral line detections, 
as we do instead. This partly explains the difference in the 
regression parameters between the two studies. Moreover, 
our measurements of [O ill] A4363 in the central 2 kpc of the 
galaxy provide additional constraints for the gradient slope 
in the inner disk. 

The linear regression to the Hll region sample, shown 
by the dashed line in Fig. 1151 is: 

12 ± log(0/H) = 8.42 (±0.03) - 0.030 (±0.008) Rkpc (2) 



Again, this is compatible with the slope of the PNe O/H 
abundance gradient. Where the PN and Hll region sam- 
ples differ is in the intrinsic scatter, which is 0.20 dex and 
0.09 dex, respectively (and easily confirmed qualitatively by 
inspecting Fig.llSfl. However, it is possible that the measure- 
ment errors for the M09 data have been under-estimated, 
because restricting the regression to our new PN sample 
the scatter is considerably smaller, 0.09 dex (with slope 
—0.027 ±0.032 dexkpc"^), confirming the suggestion made 
earlier that the abundances measured by M09 are less pre- 
cise than those we measure in our new sample. One can, 
in principle, expect a larger scatter for PNe than for Hll 
regions, due to the fact that PNe sample a significant age 
range in the evolution of galaxies. We point out that even 
in our new sample two PNe (PN 73 and PN 94) located in 
the inner 2 kpc have an appro ximately solar oxygen abun- 
dance (12 + log(0/H)Q =8.69. [Asplund et al.ll2009l ). devi- 
ating considerably (±0.2 dex) from the regression line that 
describes the abundance gradient. At the same time, PN68 
and PN 78 lie ~0.3 dex below the regression. 

We conclude that, given the uncertainties, bright PNe 
and Hll regions give a similar picture of the oxygen abun- 
dance gradient in M33. Our PNe sample provides a good 
number of auroral line-based abundances in the central re- 
gion of M33, where H 11 region abundances are more difficult 
to obtain, due to the nebular cooling. The ionized gas does 
not suggest deviations from a simple exponential gradient 
in O/H vs. radius, even in the central region of the galaxy. 

The linear reg r ession to the B-type supergiants studied 
by lUrbaneia et all (|2005h . shown in Fig. [16] by the dashed 
line, is given by: 

12 ± log(0/H) = 8.61 (±0.11) - 0.041 (±0.027) Rkpc (3) 

The O/H gradient slope from the blue supergiants is 
slightly steeper than that obtained from the ionized gas. 
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Figure 16. O/H radial gradient for PNc (red and green eircles) 
and B-type supergiants (blue star symbols). Linear regressions 
are shown by the red continuous line and the dashed blue line, 
respectively. 
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Figure 18. N/O radial gradient for PNc and Hll regions. Sym- 
bols and lines as in Fig. 1151 FN 78 was excluded from the PN 
linear regression. 
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Figure 17. N/H radial gradient for PNe and Hll regions. Sym- 
bols and lines as in Fig. 1151 PN 78 was excluded from the PN 
linear regression. 

although considering the errors the stellar gradient is still 
compatible with the gradient obtained from PNe and H 11 
regions (it is worth noting that the single star sd, R = 7.22 
kpc is largely responsible for the steeper slope). In conclu- 
sion, the oxygen abundance indicators at our disposal (PNe, 
H 11 regions and B supergiants) provide a consistent picture 
of the O/H gradient in M33. 

6.2 Nitrogen 

The radial N/H abundance gradient for both Hll regions 
and PNe is shown in Fig. 1171 The Hll regions display a 
well-defined gradient, with aslope — 0.105±0.015 dexkpc"^, 
i.e. significantly steeper than for O/H. The slope for the PNe 
is somewhat shallower (—0.060 ± 0.027 dexkpc"^), but with 
a much larger scatter (0.31 dex vs. 0.09 dex; we excluded 
the high N/H abundance PN 78 from the fit). If we interpret 
the Hll region N/H ratios as the original abundances of the 
PNe progenitors at birth, the clear offset we see between 
Hll regions and PNe in Fig. [17] is a signature of the nitrogen 
enrichment that took place during the AGB phase of the PN 
progenitors. 

The effect of nitrogen enrichment on the N/O radial 
abundance gradient is displayed in Fig. 1181 where we point 
out that the two regressions for the Hll regions and the 
PNe are virtually parallel (the slopes are —0.052 ± 0.015 



and —0.044 ± 0.025 dexkpc"^, respectively). The nitrogen 
enrichment in the PN progenitors is approximately constant, 
about 0.8 dex, across the metallicity ran ge covered by our 
samples. A similar result was found bv (Richer fc McCalll 
(2007) in the case of the SMC (where the offset between Hll 
regions and PNe was found to be ~0.7 dex), but in other 
dwarf irregulars these authors observed that the N/O ratios 
span a very wide range (~ 2 dex) for a given O/H ratio, 
so it does not seem possible to generalize our result, even 
though the N/O range in M33 appears relatively smaller. 
We also note that a few PNe have N/O ratios comparable to 
those of the Hll regions, thus their progenitors did not pro- 
duce large quantities of nitrogen, perhaps as a consequence 
of their lower masses. Although admittedly the statistics 
are rather poor, the data in Fig. 1191 where we plot N/O 
vs. O/H, suggest that this is happening preferentially at 
higher metallicity [12 -|- log(0/H) > 8.5, although this might 
depend on a poor abundance determination], and that the 
N/O ratio, and therefore the nitrogen enrichment, is larger 
at smaller O/H (a similar r esult was found, among ot hers, 
by iLeisv fc Dennefeldl Il996l and IStasinska etai\ 1 19981 : but 
note that, given the errors, the regressions in Fig. [19] are 
consistent with flat slopes). 

In Fig. [19] we also in clude the N/O r atios calculated 
for B-type supergiants by lUrbaneia et al. I (|2005l . blue star 
symbols). In these post-main sequence stars rotationally- 
induced mixing of material proces sed in the stellar core ca n 
be brought up to the photosphere (|Maeder fc Meviie3l2OO0l V 
The N/O abundance ratio, in particular, is heavily enhanced 
in B supergiants, likely by internal mixing, although addi- 
tional mechanisms (such as bi narity and magnetic fields) can 
also play an important role (|Hunter et al.l [20081 ) . B super- 
giants with even larger N/O ratios could be explained by the 
deep convection occurring during the red supergiant phase, 
altho ugh current theoret ical tracks do not extend back to the 
blue ([Hunter et al.ll2009l ). The stellar data points in Fig. [19] 
compared to the H 11 regions, show that, interestingly, the 
N/O ratio in the B supergiants is approximately the same 
as for the majority of the PNe, indicating that the nitrogen 
enrichment of present-day evolved massive stars (20-40 M0) 
is comparable to that of the progenitor stars of bright PNe 
(-2M0). 
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Figure 19. N/O for PNe, Hll regions and B-typc supergiants as 
a function of O/H. Symbols and lines as in Fig. I15l and ll6l PN 78 
was excluded from the PN linear regression. 



6.4 Radial gradient summary 

To summarize our results for this section, we report in Ta- 
ble |S] the parameters of the radial logarithmic abundance 
gradients for O, N, Ar and Ne. The number of objects used 
for the fits is included. Given the sources upon which these 
results are based, it is not surprising that we find a good 
agreement with the slopes obtained for PNe by M09 (O/H: 
-0.031±0.013 dexkpc"^; Ne/H: -0.037±0.018 dexkpc"^). 

It can be seen from Table [S] that, considering PNe 
and Hll regions separately, the slopes for the O, Ar and 
Ne abundance gradients are consistent with each other. 
In an attempt to reduce the errors for a representative 
a-element gradient slope, we have combined the abundance 
data for O, Ar and Ne, separately for PNe and Hll re- 
gions, adjusting the zero-points by the intercept amounts 
reported in Table|8] A regression to the combined data yields 




R/R25 

Figure 20. Nc/H (top) and Ar/H (bottom) radial gradients for 
PNe and Hll regions. Symbols and lines as in Fig. 1151 



6.3 Neon and argon 

The top portion of Fig. [20] displays the radial Ne/H 
abundance gradient, where we have excluded Hll regions 
with 0^'*"/(0'''^ 4- O^) < 0.5, since we have evidence that 
the ICF may be wrong for those. The slopes we obtain for 
Hll regions and PNe, -0.032 ± 0.020 and -0.023 ± 0.018 
dexkpc"^, respectively, are consistent with each other and 
with the O/H gradient slope. We must add that for the H 11 
regions the value of the slope is sensitive to the adopted 
selection in the excitation parameter: if we remove from the 
fit the Hll regions for which 0++/(0++ 0+) < 0.6, for 
example, we obtain a slope of —0.054 ± 0.018 dexkpc"^, 
more in agree ment with the resu lt obtained from Spitzer ob- 
servations bv lRubin et"al] (|2008l) . -0.058 ±0.014 dexkpc"\ 

In the bottom plot of Fig.[20]we show the radial gradient 
of the Ar/H abundance ratio for PNe and Hll regions. We 
have omitted the 4 H 11 regions with 0++/(0++ -I- 0+) >0.6. 
The gradient slope we show for the PNe is affected by the 
additional uncertainty represented by the systematic offset 
we noticed between our sample and the one by MOO. This 
uncertainty is not reflected in the gradient parameter sum- 
mary of Table |5] As the table shows, the slopes we derive 
for both types of objects agree with the values obtained for 
the oxygen abundance gradient. 



^lHgMi)=_0.033 ±0.008 dexkpc- 
dR 

for H II regions, and 

^l^i^^ = -0.018 ± 0.010 dexkpc- 
dR 

for PNe. The two slopes agree within the l-cr errors. If there 
has been an evolution of the gradient with time, it remains 
undetected with our current uncertainties. Combining the 
results, we obtain a representative slope of the radial gradi- 
ent of the a-elements in the ISM of M33: 

ISM : '^^°g("/H) ^ _o.025 ± 0.006 dexkpc"' 
dR 

We also point out that the M33 central abundances of these 
Q-elements are in good agreement between PNe and Hll 
regions. For oxygen we find a weighted central abundance 
12±log(0/H)c = 8.43 ± 0.0 2. With the adopted solar ra- 
tio [12±log(O/H)0 =8.69. lAsplund et al.|[2009t . this cor- 
responds to 0.55±0.03 times the solar value. Argon yields 
a similar result, 0.50±0.08 solar, while the neon central 
abundance equals 0.69±0.04 times the so lar value (note 
that we prefer to adopt the iLodder j |2003| solar value for 
Ar/H , which is 0.15 dex larger than the lAsplund et al.l 
I2OO9I value). The slightly discrepant result for Ne could be 
brought to a better agreement with those for O and Ar if 
we adopted a higher sola r value for Ne/0, a s advocated in 
a few recent studies (e.g. IWang fc Liull200^ ). A solar value 
(Ne/O)0 =0.22 would bring agreement between the central 
O and Ne abundances in M33, relative to solar (the most 
recent lAsplund et al.|[2009l value is (Ne/O)© = 0.17). 



7 CONCLUSION 

Our work highlights the importance of deriving auroral line- 
based chemical abundances for large samples of planetary 
nebulae and Hll regions in nearby galaxies. These observa- 
tions are essential to learn more about systematic errors in 
nebular diagnostics, the time evolution of the chemical com- 
position of star-forming galaxies, and the nucleosynthesis of 
intermediate-mass stars. Key targets for these studies are 
nearby irregular and spiral galaxies, where the spectroscopy 
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Table 8. Parameters of linear radial abundance gradients^ . 



PNe Hll regions 



Slope Correlation No. Slope Correlation No. 

Element (dex kpc~^) Intercept Coefficient objects (dex kpc~^) Intercept Coefficient objects 

O -0.013 ±0.016 8.47 ±0.07 -0.10 48 -0.030 ± 0.008 8.42 ± 0.03 -0.34 87 

N -0.060 ±0.027 8.29 ±0.12 -0.38 36 -0.105 ± 0.015 7.52 ±0.07 -0.83 21 

Ar -0.018 ±0.014 6.20 ±0.06 -0.27 36 -0.045 ± 0.016 6.27 ±0.07 -0.60 15 

Ne -0.023 ±0.018 7.75 ±0.07 -0.17 45 -0.032 ± 0.020 7.79 ±0.08 -0.27 35 



Linear relation between 12 + log(X/H) and galactocentric distance in kpc. 



of faint diagnostic lines can be carried out with the present 
generation of telescopes and instruments. 

We have obtained new deep optical spectrophotometry 
of 16 planetary nebulae in M33, mostly located in the cen- 
tral two kpc of the galaxy. For this sample we have derived 
electron temperatures and chemical abundances from the 
detection of the [O ill] A4363 line. After combining our abun- 
dances with those we obtained with th e same direct method 
from the emission line observations of iMagrini et al.l (|2009l . 
32 PNe), we have examined the behavior of nitrogen, neon, 
oxygen and argon in relation to each other in PNe and H ll 
regions, and as a function of galactocentric distances. 

The N/0 and He/H abundance ratios of most of the 
PNe of our new sample, which comprises objects in the 
brightest 1.5 magnitudes of the [O lll] A5007 luminosity dis- 
tribution, lead to a Type I classification, which, at least in 
the Milky Way, corresponds to young objects with relatively 
massive progenitors. We have found one object, PN 78, with 
an extreme nitrogen abundance, accompanied by a large he- 
lium content. 

We confirm the good correlation between Ne/H and 
O/H abundances for PNe in M33, as expected from the pri- 
mary nature of these two elements. We also find a linear cor- 
relation between argon and oxygen, but we point out that 
this is somewhat uncertain, due possibly to observational 
errors. The Ar/O ratio is consistent with the value found in 
Hll regions and in the Sun. For Ne/0, we fin d a mean value 
for PN e of 0.20, which is 18% larger than the lAsplund et all 
l|2009l ) va lue, but in agreemen t with the Ne/O vs. O/H trend 
found bv lWang fc LitJ (|2008h . These result lend support to 
the idea that at the metallicities of the PNe analyzed in M33 
these three a-elements (O, Ar and Ne) trace the abundance 
of the progenitor stars at the moment of their birth, and 
in particular that neon and oxygen have not been modified 
during the dredge-up process taking place during the AGB 
phase. 

We find no significant abundance offset between PNe 
and Hll regions at any galactocentric radius, despite the 
fact that they represent different age groups in the evo- 
lution of the galaxy. The combination of the PN and H ll 
region auroral-line abundances provide information on the 
shape of the oxygen abundance gradient in the central few 
kpc of M33. We obtain an «-element abundance gradient of 
the ISM in M33 with a slope of -0.025 ± 0.006 dexkpc~\ 
and the slopes obtained separately for H ll regions and PNe 
agree, within the current uncertainties. This is consistent 
with the conclusion by M09 that the disk of M33 has not 



experienced a significant chemical enrichment of its inter- 
stellar medium in the past few Gyr. Moreover, despite the 
large dispersion in O/H abundances for both PNe and Hll 
regions, we find no indication that towards the center of the 
galaxy PN values are systematically larger than Hll region 
ones, as would be the case if H ll region abundances were bi- 
ased towards lower values, and as could be the case if O/H 
in the central parts of M33 were really high (super-solar). 
This confirms what B supergiants already indicated for this 
galaxy. 

The growing amount of data on abundance tracers in 
M33 in the past few years has raised new questions regard- 
ing the present-day abundance gradient in this galaxy. For 
example, what is the origin of the large intrinsic dispersion 
of Hll region O/H abundances (also displayed by the PNe 
data shown here) at a given radius discovered by RS08? We 
found no correlation of the abundance residuals, as mea- 
sured with respect to the mean gradient, with a number of 
parameters, including nebular luminosity, equivalent width 
of the H/3 line, line strengths, gas density, or position in 
the galaxy disk. An analysis of additional data acquired as 
part of the M33 Metallicity Project (RS08) will hopefully 
shed some li ght on this issu e. Moreover, the study of A su- 
pergiants bv lU et al] (|2009l 'l provided metallicities near the 
galaxy center that lie above the solar value, and their de- 
rived slope, corrected to account for our adopted distance, 
is -0.084 ± 0.010 dexkpc"^ While it is true that the com- 
parison between stars and ionized gas is complicated by the 
fact that A supergiant metallicities refer mostly to the iron 
content (rather than that of the g-e lements), it is striking, 
as also pointed out bv lU et al.1 (|2009l ). that the stellar metal- 
licities are consistently near or above the solar value close 
to the center of M33, while the H ll regions samples includes 
objects with O/H ratios that extend down to 20% of the 
solar value. Future studies of the oxygen content of A and 
B supergiants in the central few kpc of M33 will help to di- 
rectly compare stellar and nebular abundances, and perhaps 
better understand the nature of the abundance dispersion. 
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